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P ■ ABSTRACT 



We present numerical simulations to describe the evolution of pre-Galactic clouds 
in a model which is motivated by cold dark matter simulations of hierarchical galaxy 
formation. We adopt a SN-induced star- format ion mechanism within a model that 
follows the evolution of chemical enrichment and energy input to the clouds by Type 
II and Type la supernovae. We utilize metallicity-dependent yields for all elements 
at all times, and include effects of finite stellar lifetimes. We derive the metallicity 
distribution functions for stars in the clouds, their age-metallicity relation, and relative 
elemental abundances for a number of alpha- and Fe-group elements. The stability 
of these clouds against destruction is discussed, and results are compared for different 
initial mass functions. We find that the dispersion of the metallicity distribution function 
observed in the outer halo is naturally reproduced by contributions from many clouds 
with different initial conditions. The scatter in metallicity as a function of age for these 
stars is very large, implying that no age-metallicity relation exists in the early stages of 
galaxy formation. Clouds with initial masses > presently observed globular clusters are 
found to survive the first 0.1 Gyr from the onset of star formation, suggesting that such 
systems may have contributed to the formation of the first stars, and could have been self- 
enriched. More massive clouds are only stable when one assumes an initial mass function 
that is not biased towards massive stars, indicating that even if the first stars were formed 
according to a top-heavy mass function, subsequent star formation was likely to have 
proceeded with a present-day mass function, or happened in an episodic manner. The 
predicted relative abundances of some alpha- and Fe-group elements show reasonable 
agreement with the observed values down to metallicities below [Fe/H] ~ —4 when the 
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iron yields are reduced relative to stellar models. The observed scatter is also reproduced 
for most elements including the observed bifurcation in [ce/Fe] for stars with low [Fe/H]. 
However, the predicted dispersion may be too large for some elements (particularly alpha 
elements) unless a limited range of progenitor masses contributing to the abundances of 
these elements is assumed. The contributions to the abundances from supernovae with 
different progenitor masses and metallicity are discussed. The results suggest that the 
low-mass end of SNell was probably absent at the very lowest metallicities, and that the 
upper mass limit for the first stars that contributed to nucleosynthesis may be <40 M Q . 
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1. Introduction 

The oldest and most metal-poor stars observed in 
the Galaxy are mainly identified as members of the 
halo population, and are believed to have formed at 
the earliest times in the history of the proto-Galaxy. 
These stars are expected to have contributed to and 
influenced the nature of all of the various popula- 
tions observed today, through their formation pro- 
cesses and the dynamical interactions that followed. 
Thus, their chemical composition and kinematic prop- 
erties can provide clues to the early epochs of Galaxy 
formation. Accordingly, many studies have been per- 
formed to identify trends in the chemical composition 
of halo stars and their connection to spatial and kine- 
matic properties (Mathews, Bazan & Cowan 1992; 
McWilliam 1995; Ryan, Norris, & Beers 1996; King 
1997; Carney et al. 1997; Preston & Sneden 2000; 
Kinman et al. 2000; Argast, et al. 2000; 2002; Car- 
retta et al. 2002; Cohen et al. 2002; 2004; Arnonc 
2004; Cayrel, et al. 2004). 

A dispersion in the abundances of heavy elements, 
varying from less than an order of magnitude for al- 
pha and iron-group elements, to more than two orders 
of magnitude for neutron-capture elements has been 
detected at very low metallicities. This has led to 
the suggestion that the lowest-metallicity stars were 
the result of mixing of the ejecta from single, or at 
most a few, supernovae, with a limited amount of gas 
in the parent clouds (Audouze & Silk 1995; Ryan et 
al. 1996; McWilliam 1997; 1998 Norris et al. 2000). 
However, the fact that the dispersion in Mg/Fe is so 
small (Arnone et al. 2004) may indicate that more 
efficient mixing may have occurred for this clement 
and perhaps other alpha elements as well (Carretta 
et al. 2002; Cayrel et al. 2004; Cohen et al. 2004). 

In addition, several studies (Lance 1988; Preston, 
Shectman, & Beers 1991; Preston, Beers, & Shectman 
1994; Majewski 1992; Kinman et al. 1994; Carney et 
al. 1996; Sommer-Larsen et al. 1997; Lec,& Carney 
1999; Chiba & Beers 2000; Chiba & Beers 2001) have 
suggested a possible duality of the halo population. 
The observed kinematics and ages of halo field stars 
suggests a hybrid formation process in which the inner 
halo may have formed by a coherent contraction in the 
manner originally suggested by Eggen, Lynden-Bcll, 
& Sandage (1962) [henceforth ELS] (see also Carney 
et al. 1996; Sommer-Larsen et al. 1997; Chiba & 
Beers 2000) while the outer halo was formed by ac- 
cretion of metal- poor fragments from nearby systems 



(e.g., Chiba & Beers 2000, and references therein). 

The most favored current theory for galaxy forma- 
tion is based on a hierarchical clustering scenario in 
which the Galaxy is assembled from cold dark mat- 
ter sub-galactic halos (e.g., Peacock 1999). These 
proto-galactic clumps are believed to have originated 
from density fluctuations in the early universe that ac- 
creted primordial gas from their surroundings. Merg- 
ing processes would lead to larger structures, and ul- 
timately, the Galaxy. Cold dark matter (CDM) mod- 
els are consistent with the observed kinematics of the 
halo of the Galaxy (Chiba & Beers 2000) and with the 
suggested hybrid formation scenario, but still require 
intensive additional observations and testing to quan- 
tify the parameters required to produce these trends 
in detail. 

Chemical evolution models that have treated these 
early stages still lack a complete picture that includes 
all chemical and dynamical effects. They are not yet 
capable of providing a comprehensive picture that re- 
produces the scatter in observed abundances at low 
metallicities. Some models in the literature (Malinie 
et al. 1993; Ishimaru & Wanajo 1999; McWilliam 
& Searle 1999; Raiteri et al. 1999; Tsujimoto et 
al. 1999; Argast et. al. 2000; 2002; Travaglio et 
al. 2000; Oey 2003) have been able to reproduce, at 
least partially, the inhomogeneities observed for el- 
emental abundances through a stochastic evolution 
scenario. These are, however, often limited for exam- 
ple to the study of only a few neutron-capture ele- 
ments, and/or metallicity independent yields, and/or 
instantaneous recycling. In the present paper we in- 
corporate metallicity dependent stellar ejecta and fi- 
nite stellar lifetimes. 

The enrichment processes that took place in these 
first structures were influenced by the initial condi- 
tions in each cloud, e.g. the initial cloud mass. While 
smaller clouds are susceptible to destruction by tidal 
forces and energetic feedback and may experience 
shorter episodes of chemical evolution, larger struc- 
tures may have continued to enrich the interstellar 
media (ISM) for longer periods of time and survived 
the merging process, contributing eventually to the 
formation of the disk. 

In this work we present a simple chemical-evolution 
model in which we evolve pre-Galactic clouds in a 
scenario that incorporates the main features of hier- 
archical galaxy formation as suggested by cold dark 
matter simulations. Our model is based upon sev- 
eral assumptions. These are: 1) In the first zero- 
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metallicity clouds, the first-generation stars are born 
as population III (Pop III), explode, and the clouds 
that survive are enriched by their expelled yields; 2) 
Star formation is induced by supernovae, mainly Type 
II (SNell); 3) Each individual supernova is taken to 
form a shell, which expands with an average velocity 
and and triggers star formation with an average effi- 
ciency e over a timescale At; 4) The shells' mass is 
the sum of the whole SNII gas and the swept-up sur- 
roundings, and the metallicity is given by a complete 
mixture of both metal contents; 5) The SN II yields 
are taken from the literature as metal dependent; and 
6) The stellar initial mass function (IMF) according 
to various conditions within the cloud. 

These assumptions, though adequate for the phe- 
nomenology of interest here, nevertheless warrant 
some caveats. For example, under the first assump- 
tion, one should in principle take care of details in 
shell fragmentation theory of shell enrichment. In the 
present work this is accomplished with effective star 
formation rates averaged over the shell. 

Although assumption 2 is an appropriate scenario 
for enhanced star formation, it neglects some hydro- 
dynamic negative self-regulation. For example, the 
hydrodynamic expansion of a SN bubble depends on 
the external density. In the hydrodynamic Sedov 
phase (Sedov 1946; Ostriker & McKee 1988; Chioffi, 
McKee & Bertschinger 1988) energy is lost due to ex- 
pansion because the expansion velocity decreases with 
time, and also depends on density as p" 1 / 5 . When 
the ambient density is larger at an earlier times, both 
radius and velocity are less, so the sweep-up of sur- 
rounding gas could be smaller than assumed here. 
Furthermore, if subsequent stars are formed in as- 
sociations, SNII bubbles of massive stars would be 
expected to overlap and to form a so-called super 
bubble. This means that the shells' total surface is 
smaller, and the amount of swept-up gas is reduced. 

In our model all of these details are absorbed into 
simple average properties of the clouds and shells. 
Nevertheless, these assumptions are adequate to de- 
scribe the bulk of the observations, as we shall see. 
In particular, we incorporate the stochastic nature 
by which the initial conditions including total mass, 
and metallicity must have been distributed among the 
different clouds according to a hierarchical clustering 
scenario. This random nature is further extended dur- 
ing the early stages of evolution by the dependence of 
the chemical yields produced in SN events on the ini- 
tial mass and metallicity of their progenitors. There- 



fore, we consider clouds of different initial masses 
and form stars according to four different IMFs. We 
also consider metallicity- dependent chemical yields of 
a number of iron group and alpha elements allowing 
for the products of SN events to change with time. 

Assumption 4 is another very crucial assumption, 
since it strongly enters the metallicity results. The 
interstellar medium of galaxies includes a hot metal- 
rich gas phase that also carries these elements and 
could evaporate them from the clouds. Hence, this as- 
sumption is relevant for element abundances and the 
timescale of their enrichment. Indeed, dwarf galax- 
ies may require inefficient mixing of supernova ejecta 
with the ISM to achieve observed low oxygen abun- 
dances and other low effective yields. We note, how- 
ever, that although the full yields of SNell may not 
be instantaneously incorporated and homogeneously 
mixed into the cool gas, in our model this loss of SN 
mixing efficiency can be compensated for by an in- 
crease in the star formation efficiency. Hence, for our 
purposes an efficient mixing of the supernova ejecta 
is an adequate working hypothesis. We return to this 
point below. 

The specific form for SN-induced star formation in 
our model is based upon that proposed by Tsujimoto 
ct. al (1999). That is, star formation is triggered by 
SN events, and the material from which the next gen- 
eration of stars forms is the result of complete mixing 
of the SN ejecta with the gas swept up from the ISM 
by the explosion. This produces low-metallicity stars 
with elemental abundances resembling the ejecta of 
high-mass SN progenitors, as suggested by the data. 

We show that this simple stochastic model is ca- 
pable of reproducing the observed scatter in abun- 
dances, relative to iron, at low metallicity, as well as 
the observed shift in trends at higher metallicities af- 
ter the chemical products of stars mix with the ISM. 
The metallicity-dependent chemical yields of Woosley 
& Weaver (1995, hereafter WW95) used in this model, 
have been used previously in detail by Timmes et al. 
(1995), who showed general agreement with the ob- 
served relative abundances in the Solar neighborhood 
for most elements. They did not, however, attempt 
to apply these chemical yields to a stochastic, CDM- 
motivated model for the halo, and they only did their 
comparison for [Fe/H] > —3. 

In the present work we test these yields with our 
model for several alpha and several iron-group ele- 
ments down to metallicities below [Fe/H] ~ —4.0, 
near the present observational limit of metallicity in 
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the Galaxy. 

2. The Model 

The main purpose of the present study is to ex- 
plore the viability of a simple schematic hierarchi- 
cal clustering paradigm for galaxy formation. In this 
picture, primordial fluctuations lead to the gravita- 
tional collapse of subsystems forming sub-galactic ha- 
los comprised of both baryons and dark matter. Such 
subsystems then merge and accrete diffuse matter to 
form proto-galactic halos. We follow the idealized 
chemical evolution of the baryon component of in- 
dividual clouds, assuming all stars formed initially in 
proto-Galactic clouds via SN-induced star formation 
(e.g., Tsujimoto et al. 1999). We utilize metallicity- 
dependent chemical yields and follow the evolution of 
several alpha- and Fe-group elements. These are com- 
pared to the observed abundances of individual stars 
in the halo and thick disk of the Galaxy. 

Our model consists of evolving individual clouds 
in the mass range of 10 5 to 10 8 M©. We follow the 
chemical evolution of each cloud as a one-zone closed 
box (Tinsley 1980). The evolution of gas and stars 
is then explicitly evolved in time. We assume that 
all halo field stars were formed originally in these 
proto-Galactic clouds and were later dispersed by 
tidal forces. By varying the total masses and the star 
formation histories of these clouds, we attempt to ac- 
count for the trends and scatter in elemental abun- 
dances observed in metal-poor stars of the halo and 
thick-disk populations. 

Simulations show that as the baryonic component 
of these clouds cools (mostly via molecular hydrogen) 
they gravitate to the center of the dark-matter poten- 
tial to form central cores (Abel et al. 1999; Norman 
et al. 2000). These cores are believed to form stars 
at their centers, the exact nature of which is still un- 
certain (Silk 1983; Carr et al. 1984). We assume that 
each cloud starts with primordial material, and ulti- 
mately forms one massive star at its center. This Pop- 
III star later explodes as a SN event, triggering the 
formation of higher mctallicity (population II) stars 
in its high-density shell. Subsequent star formation 
progresses in the shells of later SN events. 

Star formation has been linked to the dynamical 
evolution of molecular clouds in many studies. Obser- 
vations of nearby molecular clouds suggest that star 
formation is triggered by expanding shells driven by 
stellar winds and repeated SN events in galaxies (Ya- 



maguchi et al. 2001; Hartmann 2002). When these 
shells become unstable and fragment, they may form 
molecular clouds and stars (Ehlcrova et al. 1997; 
Blondin et al. 1998). The conditions under which 
star formation is triggered due to the gravitational 
collapse of expanding shells has been investigated by 
Elmegreen et al. (2002) with a 3-D numerical simu- 
lation of the conditions in the ISM as the expanding 
shell injects its energy. 

In all of these studies, stars are assumed to form 
out of the gas in the ISM and reflect its metal con- 
tent at the time of their formation. Although this 
treatment of the star formation process in galactic 
chemical evolution models has been very successful in 
reproducing the general features of the chemical com- 
positions of stars and HII regions in the solar neigh- 
borhood (Matteucci & Greggio 1986; Yoshii, Tsuji- 
moto, & Nomoto 1996), it cannot be applied to the 
early stages of the galaxys evolution (Shigeyama & 
Tsujimoto 1998), as it fails to reproduce the observed 
abundance patterns in extremely metal-deficient stars 
(Mcwilliam et al. 1995; Ryan et al. 1996). Therefore, 
an alternative scenario was proposed by Tsujimoto et 
al. (1999) in which stars are born in the dense shells 
formed by the sweeping up of ISM during a SN explo- 
sion. Ryu & Vishniac (1987) suggest that such dense 
shells which form behind the radiative shock front 
are dynamically overstable and are broken into a few 
thousand fragments in which self gravity is unimpor- 
tant. Nakano (1998) suggests that some of these frag- 
ments will dynamically contract and eventually form 
stars. Such a SN-induced star formation scenario has 
been shown to fit the observed relative abundances at 
early times better than models where SN ejecta mix 
completely with all the gas in the ISM (Nakamura 
et al. 1999). We adopt this scenario with a new set 
of equations consistent with our use of metallicity- 
dependent yields. The metallicity of each shell is cal- 
culated at the time of its formation and stored for use 
in subsequent generations. 

This scenario leads to a difference between the 
abundances observed in stars and those of the ISM 
at the time they were formed. This differs from most 
Galactic chemical evolution models, in which the stel- 
lar and gas abundances are constrained to be identi- 
cal. The material from which the next generation of 
stars forms, following a given SN explosion, is the re- 
sult of the complete mixing of the SN ejecta with the 
gas swept up from the ISM by the explosion. This 
star formation scenario will produce low-metallicity 
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stars with elemental abundances resembling that of 
a single high-mass supernova. Our model follows the 
evolution of gas in the clouds and the evolution of 
individual elements in the ISM. We also follow the 
chemical evolution of stars and derive the anticipated 
present-day metallicity distribution function for halo 
stars. 

2.1. The Calculation 

We performed numerical solutions of the chemical- 
evolution equations for each cloud, deducing the evo- 
lution of gas mass, M g , and mass of individual el- 
ements, Mj, in finite steps of time. At each time 
step, the equations are integrated over a grid of stellar 
masses, ranging from 0.09 to 40 M Q . In the context 
of this model, we desire a time step, At, longer than 
the time required for a SN shell to form and disperse, 
but short enough to be sensitive to the short lifetimes 
of massive stars. For example, the expansion time 
of a SN event is of the order of 0.01 Myr (Binney & 
Trcmainc 1987), while the lifetime of the most mas- 
sive stars included in our model is ~ 6 Myrs. Since the 
sensitivity to the lifetimes is very crucial at the earli- 
est times, the time step was taken to be ~ 1 Myr up 
to the point when the least-massive SNe would have 
contributed (~20 Myrs). Thereafter, At was given a 
larger value of 25 Myrs. 

We take into consideration the main sequence life- 
times of these stars, relaxing the instantaneous recy- 
cling approximation. This is of fundamental impor- 
tance when dealing with timescales as short as the 
halo formation time, and when treating intermediate- 
mass elements, which receive significant contributions 
from progenitor stars with lifetimes comparable to the 
halo formation time. We adopt the mass-dependent 
lifetimes used by Timmes et al. (1995), i.e. the life- 
times of Schaller ct al. (1992) were used for stars less 
massive than 11 M@, while the lifetimes given by the 
stellar evolution calculations of Weaver & Woosley 
(1994) were adopted for more massive stars 

For a given cloud, when the first shell forms at 
t = 0, the rate at which mass goes into stars in the 
cloud at that time is 

ij(t = 0) = eM sh (m,0)/At , (1) 

where e is the star formation efficiency which gives 
the mass fraction of the shell that is used up in form- 
ing stars. It is treated as a free parameter in our 
model. Changing this parameter can only have an 



effect on the time scales, but will not alter the over- 
all trends calculated for the elemental abundances or 
the MDFs. The e parameter was adjusted in this 
model to reproduce the shift in slope of the elemental 
abundance ratios for three iron-group elements below 

[Fe/H] 2.4 observed by McWilliam et al. (1995). 

It was also made to be consistent with the obser- 
vation that the age-metallicity distribution changes 
for extremely metal- poor stars with [Fe/H] < —2.4. 
M s h(m,t) is the mass of the shell formed in a SN ex- 
plosion at time t, with the progenitor mass being to, 
and is given by : 

M sh (m,t) = E j (m,Z) + M sw , (2) 

where M sw is the mass of the interstellar gas swept 
up by the expansion. Since the explosion energy of a 
core-collapse supernova depends only weakly on the 
progenitor mass (WW95; Thielemann et al. 1996), 
this quantity is taken to be constant with a value of 
5x 10 4 M (Ryan et al. 1996; Shigeyama & Tsujimoto 
1998; Tsujimoto et al. 1999). Ej(m, Z) is the mass of 
all the ejected material from the SN with a progenitor 
mass to and metallicity Z. [Note that if not all of the 
supernova ejecta is mixed into the ISM, then this can 
be compensated in out model by the star-formation 
parameter e.] For the first SN event the metallicity Z 
is zero for a population III star, and to is toi. For later 
generations, the metallicity of a star is calculated from 
the metallicity of the shell from which it was formed, 
Z s h(m,t). This is calculated as a function of time, 
and is given by: 

Z sh (m,t) = l-[x» h (m,t)+x? h e (m,t)} . (3) 

Here, x^ h and x^ are the fractions of H and He in 
the shell, respectively, and are given by: 

x^ h (m,t) = [y H (m,t)+xf as (t)M sw }/M sh {m,t) , (4) 



x? h e (m, t) = [y He (m, t) + xf a e s (t)M sw ]/M sh (m, t) , 

(5) 

where y H and y He are the mass of ejected H and 
He, respectively, from the explosion. We will refer to 
these quantities as the "yields" . They are metallicity 
dependent and therefore are functions of time. The 
quantities x¥ as and x^ s are the fractions of H and 
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He in the interstellar gas at the time of formation of 
the shell t, respectively. This fraction is given for any 
element x gas by: 

x i gas (t)=M i (t)/M g (t) . (6) 

For the first explosion, yi(mi,t\) is the yield of 
elements i from stars of mass m\ and metallicity zero. 
The subsequent star formation rate (SFR), ij){t), is 
derived from a sum over all shells that form at time t. 
It thus depends on the SFR at the time the progenitor 
star formed, tp(t — T m ), where r m is the lifetime of the 
progenitor. 



4>(m) 



m 



ip(t>0) = / eM sh ( mi ,0) 

J max(mt.lO) 

x ip(t — T m )dm 
where <j){m) is the IMF, normalized such that 



/' 

J m 



4>(m)dr 



1 . 



(7) 



(8) 



The quantities m u and mi denote, respectively, the 
upper and lower mass limits for the IMF. The lower 
mass limit for stars that produce SN events of type-II 
will be taken as 10 M©. The quantity mt is the mass 
of a star with lifetime equal to t, measured from the 
time of formation of the first shell. Since the ejected 
mass from a SN changes with the metallicity of the 
progenitor, the quantity Ej(m, Z s h) in Eq. (2), must 
be replaced by an average over all stars of mass m 
and different values of Z s h , which explode at a given 
time t: 



(M ej ( 



.*)>= / 

J max{r 



dm' 



m(t-T m ),W} 

{4>{m')/m')m ej {m,Z sh {m',t- r m ) 



, (9) 



Here, m' is the mass of the progenitor which produces 
the shell from which m is formed, and m e j (m, Z s h (m', t- 
T m )) is the mass ejected from the star with mass m 
and metallicity Z s h(m! , t— m ). The quantity (t — r m ) 
denotes both the time of formation of a star of mass 
m and the death of the star m' . 

The same applies for the yields of individual ele- 
ments, yi 



(yi{m,t)) 



J max{ 



dm'i 



m(t-r m ), 10} 



(f){m')/m' 



x y? (m,Z sh {m',t-T m )) . (10) 

and y\ 3 (m, Z s h(m' , t — r m )) is the mass of the element 
i ejected from the star with progenitor mass m and 
metallicity Z s h(m', t — r m ). 

The change in gas mass with time is then given by: 



dM g 
dt 



-m + / 



raax{va t ,'mi) 



dm[(j)(m)/m] 



xM ej (m,t)tp(t-T m ) . 



(11) 



The first term in this equation is the SFR, equal to 
the amount of gas going into forming stars at time t. 
It is derived from Eq. (7) . The second term accounts 
for the enrichment process by all stars whose life ends 
at time t, including the whole possible range from m t 
to m u . The change in the mass of element i in the 
gas is given by: 



dMi 
~df 



L 



max(m tl W) 



dm[(p(m)/m] 



xeM sh (m,t)xi(m, t)tp(t - r m ) 
dm[(j)(m)/m] 



max(m tl mi) 

xyi(m,t)il>(t-T m ) . 



(12) 



The first term in this equation gives the rate at which 
element i is incorporated into stars at time t, whereas 
the second term represents the rate at which it is 
being added to the ISM by stars. By substituting 
Eqs. (7)-(9) into equations (10) and (11), this set 
of intcgro-diffcrential equations is solved numerically. 
An auto-regressive computation of the star formation 
rate and the metallicity of each shell that forms was 
employed, and previous values were substituted in 
equations (7) - (11) at each time step. At every time 
step, the metallicities of all shells formed are calcu- 
lated and stored for use at later times. Each shell 
is recorded and the age-metallicity relation is con- 
structed. The evolution of the metal content in both 
the gas and and the stars are followed as a function 
of time. 

2.2. Metallicity Dependent Chemical Yields 

2.2.1. SNII Yields 

In this model we have utilized the stellar yields 
calculated by WW95 for high- mass stars. These 
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are more-or-less consistent with the yields of Porti- 
nari, Chiosi & Brcssan (1998). Our study is there- 
fore complementary to the models of Argast et al. 
(2000; 2002) who made a similar stochastic study, 
but relied upon the yields of Thielemann, Nomoto, & 
Hashimoto (1996) and Nomoto et al. (1997). We use 
finely spaced mass and metallicity grids in our simu- 
lation, and linearly interpolate the yields in both di- 
mensions. The yields of every SN event are mass- and 
metallicity-dcpcndcnt, and the metallicities are calcu- 
lated for all stars at the time of their formation and 
stored for use at later times. WW95 included stars 
with progenitor masses ranging from 11 to 40 M©, 
and metallicities from zero to solar, and to calculate 
yields for elements from H to Zn. This permits better 
limits to be set on Galactic chemical evolution models, 
and enables comparisons between the results of sim- 
ulations and the empirical data from high-resolution 
spectral measurements in stars. 

The uncertainties associated with the use of such 
supernova yields have been critically summarized in 
Argast et al. (2000; 2002). In particular, the syn- 
thesized yields can be sensitive to various aspects of 
both the stellar nucleosynthesis models and the super- 
nova explosion mechanism and energy released. As a 
result one needs to be cautious in inferring quanti- 
tative conclusions regarding detailed elemental abun- 
dances. Nevertheless, qualitative trends in elemen- 
tal ratios can be studied. By exploring the trends 
based upon independent models utilizing different in- 
dependently derived yields it is at least hoped that 
the present study will help to clarify the underlying 
issues in early Galactic chemical evolution. In this 
sense it is important to compare the present study 
with that of Argast et al. (2000). 

It is useful to consider some of the differences be- 
tween our yields and those employed by Argast et 
al. (2000). While the yields presented by WW95 
cover stars up to 40 M Q , other studies, such as Thiele- 
mann, Nomoto & Hashimoto (1996) obtain yields for 
stars up to 70 M Q . Although the fate of stars more 
massive than 40 M© is still uncertain, recent studies 
have suggested that metal-free stars with masses be- 
tween 35-100 M collapse into black holes, while stars 
with masses between 10-35 M can explode as SNe- 
II, and contribute to the enrichment process ( Heger, 
Woosley, & Waters 2000; Heger et al. 2002). In addi- 
tion, including such massive stars will have very little 
effect on the results of a chemical-evolution calcula- 
tion when adopting a Salpeter IMF, since very few 



stars will form that are more massive than 40 Mq. 
The metallicity dependence of the yields of some ele- 
ments, such as oxygen, is very strong for such massive 
stars. This dependence will influence greatly their 
evolution at early times. Since oxygen is not included 
in our study, we believe that the inclusion of these 
massive first stars is not crucial here. Therefore, the 
stellar mass range encompassed by the calculations of 
WW95 is reasonable. This is also consistent with the 
findings of Samland (1998). 

In contrast to the stellar yields calculated by Thiele- 
mann et al. (1996), in which they use a constant solar 
progenitor metallicity, the Woosley & Weaver yields 
have the advantage of covering the whole metallicity 
range from zero to solar for the progenitor stars. Al- 
though their predicted yields of heavy elements are 
not affected greatly by small changes in the initial 
metallicity of the progenitor (for Z/Zq > 0), there 
are large differences between the predicted yields for 
stars with metallicity Z/Zq > and those for stars 
with Z/Zq = 0. This has a significant effect on the 
metallicity distribution of extremely metal-poor stars, 
and on the stars of the next generation that will form 
in their shells. 

WW95 also emphasized that, in addition to the ini- 
tial mass and metallicity, the energy of the explosion 
is an important parameter in determining the yields 
from a SN event, The energy of the explosion deter- 
mines the mass cutoff between the part of the synthe- 
sized elements that will be ejected and elements in the 
deeper layers that will fall down onto the core after 
the explosion. They distinguished different models A, 
B, and C, for stars more massive than 25 M©. These 
correspond to three different values for the initial ki- 
netic energy of the "piston" (a theoretical construct 
used to simulate a real explosion). The larger this 
energy, the more heavy elements escape from the ex- 
plosion. We adopt model B, an intermediate value for 
the energy. 

2.2.2. Intermediate- Mass Yields 

For intermediate stellar masses we adopt the stan- 
dard nucleosynthesis yields of Renzini & Voli (1981). 
Although other tracks are available (e.g. Portinari 
et al. 1998), the Renzini and Voli (1981) tracks 
are sufficient for the present study, which is almost 
unaffected by the yields of intermediate mass stars 
other than through the overall evolution of metallic- 
ity. The Renzini & Voli (1981) yields are also the 
same as those adopted in earlier works (e.g. Timmes 
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et al. 1995) and so are useful for comparison of the 
present stochastic model with the results of that ear- 
lier continuous-star-formation model. 

The Renzini & Voli (1981) models included stars 
in the mass range of 1 < m < 8M Q . They calculated 
the yields for two different metallicities, Z = 0.004, 
and Z = 0.02. In the current model, we interpolated 
between these two values and used Z = .004 yields for 
lower metallicities. Intermediate-mass stars mainly 
eject H, He, C, and N into the ISM at the end of their 
lives. Thus they mainly only affect the metal/H ratio. 

2.2.3. SNIa Yields 

Intermediate-mass stars mostly end their lives as 
carbon-oxygen rich white dwarfs. If they are members 
of a binary system these dwarfs may accrete mass 
from the remaining member to the point where their 
mass exceeds the Chandrasekhar limit. At this point, 
the star becomes unstable, causing a thermonuclear 
Type-la SN event. For these events we use the yields 
calculated by Thielemann et al. (1986). 

The appearance of the Fe-rich material contributed 
by Type-la SN is delayed by the time required by a 
white dwarf to evolve to an explosion. This delay 
time has been estimated to be of the order of 1 Gyr 
(Greggio & Renzini 1983; Matteucci & Greggio 1986; 
Smecker & Wyse 1991; Mathews et al. 1992; Yungel- 
son & Livio 1998; Tantalo & Chiosi 2002). There- 
fore, these events will only have a significant effect on 
clouds that evolve and self-enrich on a timescale of 
more than 1 Gyr. 

Disrupted clusters disperse and permeate the halo. 
The dispersed remnants can produce isolated SNela 
but should not significantly affect the cloud abun- 
dances of interest in the present study. Only SNIa 
events which occur within the clouds affect the present 
study, which can be the case for clouds that are suf- 
ficiently massive to retain their gas. Less massive 
structures are expected to be destroyed due to the en- 
ergy input from SN of Type-II events during the first 
Gyr of the onset of star formation. Stars less massive 
than 1 Mq serve in this model only as reservoirs of 
gas mass, since they do not evolve significantly during 
the considered time. 

2.3. The IMF 

In addition to the type of mixing processes follow- 
ing the explosion, the enrichment of the halo depends 
on the type and number of SN events that took place. 



Therefore, the central role that the IMF plays at these 
early times cannot be ignored. For a SN-induced star 
formation model, the SFR depends on the rate of 
death of massive stars at a given time. In other words, 
the SFR today depends on the SFR at previous times. 
Therefore, the IMF of the first stars (perhaps more 
properly referred to as the First Mass Function, FMF) 
will affect greatly the subsequent evolution of these 
clouds. 

In this regard it is of interest that several re- 
cent studies have suggested that Pop-Ill stars where 
probably very massive (Elmegreen 2000a; Bromm et 
al. 2001b; Bromm et al. 2002; Abel et al. 2002; 
Kroupa 2002; Nakamura & Umemura 2002), or may 
have had a very massive component (Nakamura & 
Umemura 2002) compared to modern stellar popula- 
tions. 

Although the shape of the IMF that dominated 
in primordial conditions is still uncertain, there is 
compelling circumstantial evidence pointing toward 
an early IMF biased towards more massive stars. 
This evidence includes the paucity of metal-poor stars 
in the Solar neighborhood [i.e. the G-dwarf problem 
(Schmidt 1963; Larson 1998)], observations indicat- 
ing a rapid early enrichment of the ISM by heavy el- 
ements produced in the explosions of high- mass stars 
(Mushotzky et al. 1996; Loewenstein & Mushotzky 
1996), followed by slower enrichment at later times, 
represented by a flat age-metallicity relation in the so- 
lar neighborhood (Larson 1986). Further evidence is 
found in the high abundances of heavy elements in the 
hot gas trapped in the potential wells of rich clusters 
of galaxies. The heavy-element abundance is larger 
than what would be predicted by a standard present- 
day IMF (Zepf & Silk 1996 and references therein; 
Larson 1998). 

Moreover, several studies of the central regions of 
star burst galaxies have suggested that the IMF is 
strongly biased towards high-mass stars (Doane & 
Mathews 1993; Rieke et al. 1993; Doyon, Joseph, & 
Wright 1994). Observations of nearby ellipticals show 
that [Mg/Fe] increases above the solar value with in- 
creasing galaxy luminosity (Worthey, Faber & Gon- 
zalez 1992; Davies, Sadler , & Peletier 1993) , and of 
distant galaxies with unusually red colors or strong 
4000 A breaks (Chariot et al. 1993). These results 
are suggestive of an IMF that was once weighted to- 
ward massive stars, perhaps during the early merger 
that was the hallmark of a forming elliptical. 

A top-heavy IMF during the initial stage of ellipti- 
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cal galaxy formation has also been proposed as a way 
to account for the large iron masses in galaxy clus- 
ters (Elbaz, Arnaud, & Vangioni-Flam 1995). Recent 
observations of the elemental abundances of the intr- 
acluster gas point to Type II SN as the source of en- 
richment, providing strong support for models with 
a top-heavy IMF during the formation of elliptical 
galaxies (Loewenstein & Mushotzky 1996). 

In addition, the elemental abundances observed in 
halo stars suggest that they were made out of gas 
enriched only by SNell (Nissen et al. 1994; Ryan et 
al. 1996; Tsujimoto et al. 1999; Norris et al. 2000). 
Finally, recent theoretical calculations (Nakamura & 
Umemura 2002) of star formation in extremely metal- 
deficient clouds suggests that, in such conditions, the 
mass function is peaked at the high mass end (10-100 
M Q ) and is deficient in sub-solar mass stars. 

The Salpetcr power-law IMF represents the stellar 
distribution in the solar neighborhood very well down 
to 1 M Q . At lower masses the form of the function 
is not clear due to the difficulty of obtaining a mass- 
luminosity relation for such faint stars. However, it is 
believed that it must decline rapidly and flatten below 
0.1 M Q in order to explain the paucity of observed 
brown dwarfs compared to the predictions (Basri & 
Marcy 1997; Larson 1998; Elmegreen 2000b; Kroupa 
2002). 

Empirical studies of the IMF have been conducted 
recently in different environments including clusters 
and associations in the Galaxy and the Magellanic 
Clouds (von Hippel et al. 1996; Hunter et al. 1997; 
Hillenbrand 1997; Massey & Hunter 1998). Although 
all of these studies support a Salpeter IMF with a 
slope in the neighborhood of ~ 1.35 for stars more 
massive than a solar mass, some different slopes have 
been suggested (Scalo 1998; Massey & Hunter 1998). 
A top-heavy IMF for the first stars has also been sug- 
gested by a number of authors (e.g., Larson 1998; 
Bromm et al. 2001a). Two functional forms were 
proposed by Larson (1998). 



dN 
dlogm 



a (1 + m/mi) 



-1.35 



and 



dN 



oc to cxp{— mi/m} 



(13) 



(14) 



dlogm 

These both approach the power law with a Salpeter 
slope at high masses and fall off at the low-mass end. 

Eq. 13 falls off asymptotically to a slope zero at 
the low end. Eq. 14 has a peak at to = m;/1.35, 



and falls off exponentially with increasing negative 
power at lower masses. The characteristic mass scale, 
to;, has a value of 0.35 M Q for the present-day so- 
lar neighborhood. This mass function was used by 
Hernandez & Ferrara (2001), who suggest that the 
IMF for population III stars was strongly weighted 
towards high masses at redshifts 6 < z < 9. They 
infer a value for mi of 12.2 at redshift z ~ 9, and 
metallicity [Fe/H] < — 2.5, based on number counts of 
metal-poor stars from the HK survey (Beers, Preston, 
& Shectman 1985; 1992). Both functions were tested 
in this model and compared to the simple power law. 
The mass scale was chosen to be a linear time depen- 
dent quantity of the form to; = 13.2 — 0.9t. The range 
of the IMF is taken from 0.09 to 40 M Q . 

Figure 1 shows the different mass functions ex- 
plored in this study. The first is the simple Salpeter 
function with index = 1.35. The second is a more 
moderate high-mass biased Salpeter-like function that 
shows a slower decline towards high masses than the 
simple function, but still steeper than the top-heavy 
functions of Larson. The two functions from Larson 
behave differently at low masses. While the first pro- 
duces a fair amount of low-mass stars and declines 
slowly at the massive end, while the second function 
produces very few low-mass stars and peaks at an in- 
termediate value before it declines exponentially 

2.4. Initial Conditions 

The baryonic masses of globular clusters observed 
in the Galaxy range between 10 4 and 10 6 M Q , while 
the minimum total mass of dwarf galaxies in the Lo- 
cal Group arc - 2 x 10 7 M Q (Mateo 2000). Therefore, 
we choose initial total cloud masses in the range (10 5 - 
10 8 ) M Q . In this mass range, the clouds are sensitive 
to stellar feedback, both in terms of chemical enrich- 
ment or energetics. We start with gas in a single 
phase and uniform density, and choose the volume of 
each cloud such that we maintain a reasonable initial 
density in the clouds (~ 0.25 particles cm" 3 ) corre- 
sponding to a baryon surface density of about 0.01 
M pc~ 2 . This condition is necessary in order to cal- 
culate the initial potential energy of the cloud, and 
its lifetime prior to destruction by SN energy input. 

The dark matter component of the clouds is taken 
to be ~ 3 times the baryon mass, consistent with ob- 
servations of rich galactic clusters (e.g., Mushotzky 
et al. 1996). The baryonic component of the clouds 
are taken to be composed initially of 77% (by mass) 
Hydrogen, 23% Helium, and a zero metallicity. The 



10 



first massive star forms at the core and initiates sub- 
sequent events of SN-induced star formation. This 
method of star formation produces local inhomo- 
geneities in the clouds as soon as the first SN event 
takes place. Each SN event will produce a group of 
stars reflecting the abundance pattern of that partic- 
ular SN progenitor, since the stellar yields of a SN are 
different for different progenitor masses. 

It has been shown (Audouze & Silk 1995) that mix- 
ing timescales in the early halo are sufficiently long 
that chemical inhomogeneities in the gas would not 
be erased on the timescale over which early genera- 
tions of stars form. Therefore, we do not allow mixing 
between shells. Still, we assume instantaneous mix- 
ing of the shell material left behind after the explosion 
with the gas in the ISM. 

Our clouds are self-enriched over a timescale of ~ 1 
Gyr. The timescale for the formation of the early halo 
is estimated to be of the order of a few Gyr (Bekki 
& Chiba 2000), rather than 10 8 yrs as suggested by 
ELS, although most of the processes that contributed 
to the formation of the earliest generations of stars 
will take place in the first 0.5 Gyr after the initiation 
of the first enrichment events. 

3. Results and Discussion 
3.1. Stability of the Clouds 

Ignoring any tidal effects, heating and cooling, and 
any inhomogeneities in the gas, and taking the energy 
input per SN event Esn to be 10 51 erg, we find that 
under the adopted initial conditions, clouds with total 
mass <3 x 10 6 M Q (baryon mass <7 x 10 5 M Q ) do 
not survive the first SN explosion, since their total 
gravitational potential energy is less than Esn- For 
more massive clouds, the survival time depends on 
the SFR chosen and on the form of the IMF. 

The initial density has a large effect on the stability 
of a cloud, and we must keep in mind that star- form- 
ing regions become cooler and denser as time passes. 
As they form stars, the clouds are heated by the ther- 
mal input from massive stars. Therefore, such a sim- 
ple treatment can only give a rough estimate of the 
fate of these clouds. Table 1 shows the survival times 
of our model clouds against destruction, as a function 
of initial mass, for three different mass functions. 

It is clear from the table that the fate of clouds 
which survive the first explosion varies greatly de- 
pending on the type of the IMF chosen. The Salpeter 



function produces very few massive stars, which al- 
lows enrichment to continue in all clouds more mas- 
sive than 5 x 10 6 M Q for times longer than 10 Gyr. 
At the same time, clouds as massive as 10 8 M Q do 
not survive more than ~ 0.1 Gyr with a top-heavy 
IMF. This leads to the suggestion that even if the 
first stars were formed according to a top-heavy IMF, 
subsequent star formation must have involved an IMF 
similar to the present day mass function. 

This result is in agreement with recent studies con- 
cerning the type of IMF expected to govern mass for- 
mation in the past and present (Larson 1998; Her- 
nandez & Ferrara 2001; Nakamura & Umemura 2001; 
Nakamura & Umemura 2002). Larson (1998) sug- 
gested the existence of a mass scale in the star- 
formation process which varies with cloud conditions 
such as pressure and temperature, and therefore is 
time dependent. If the temperature in star- forming 
clouds was higher at early times, this will have the 
effect of increasing the mass scale and the relative 
number of low-mass stars formed at early times will 
decrease. He also suggested another possible form of 
IMF with a universal power-law form at large masses, 
but that departs from this power law below a charac- 
teristic time-dependent mass scale. 

Empirical evidence has supported this picture, with 
a power-law form of the IMF down to one Solar mass 
(see §2.3). Recent chemical evolution models have 
tested the time varying IMF in our galaxy as a so- 
lution to the G-dwarf problem (e.g., Martinelli & 
Matteucci 2000). These authors have found that an 
IMF with two slopes, and a time- dependent shape 
at the low mass end, is required to reproduce con- 
straints other than the G-dwarf problem. They also 
show that such a function is still not able to repro- 
duce the properties of the Galactic disk and suggested 
the inclusion of radial flows to reproduce the observed 
trends. Invoking different slopes of the IMF in order 
to explain the G-dwarf problem is a long-standing way 
to address this problem which has fallen out of favor 
(see above). More recently, however, dynamical ef- 
fects, both of the infall of low-metallicity gas, and 
radial outflow of hot metal-rich gas (e.g., Samland et 
al. 1997) have provided a somewhat more plausible 
mechanism to produce these effects. Nevertheless, the 
question still remains as to whether the mass function 
has varied with time and needs to be addressed fur- 
ther by including dynamical effects in the models of 
galactic evolution. For the purposes of the present 
study a simple time-dependent IMF is adequate. 
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Table 1 
Survival Times of the Clouds 



M c i oud (M ) TO- 1 ' 35 (1 - e -W2) m -1.35 m -1.35 e -m,/4 



3 x 10 6 > 10 Gyr 

5 x 10 6 > 10 Gyr 9 Myr 8 Myr 

1 x 10 7 > 10 Gyr 80 Myr 12 Myr 

1 x 10 8 > 10 Gry 200 Myr 90 Myr 



Our model is consistent with the clouds being the 
environment out of which the globular clusters first 
formed. It is necessary for the efficiency e of star 
formation be very high for globular cluster forma- 
tion. The gravitational binding of the star cluster 
that is formed must therefore dominate, because gas 
expulsion did not lead to the cluster's disruption. The 
cloud mass which survives disruption (~ 7 x 10 5 M ) 
is consistent with observed globular cluster masses, 
particularly if some evaporation has occurred to the 
present time. In our model, sufficiently massive GCs 
could have been self enriched. Indeed, a recent more 
detailed model by Parmcnticr (2004) for the forma- 
tion of globular-cluster systems suggests that they 
were indeed self enriched. A caveat, however is that 
if the efficiency ijof star formation is very high, as 
e.g. for globular cluster formation, the gravitational 
binding of the star cluster that is formed must domi- 
nate because gas expulsion did not led to the cluster's 
disruption. 

Alternatively, globular clusters could have been 
formed initially in more massive stable systems such 
as dwarf galaxies. The homogeneous chemical com- 
position of stars within most globular clusters seems 
to suggest that they formed out of gas that was al- 
ready pre-enriched and well mixed inside structures 
that were more stable against destruction. The glob- 
ular cluster population of the halo might then be ex- 
plained as the result of accretion events during tidal 
interactions with other galaxies. This was suggested 
previously by Armandroff (1993) and Larson (1999). 

As mentioned previously, our calculation of the 
lowest mass of a cloud that can survive the first SN 
event is very simplistic and should, perhaps, not be 
taken too seriously on this issue of globular cluster 



formation and survival. Among other things, our 
model ignores hydrodynamical effects and inhomo- 
geneitics/fragmentation in the gas. More precise un- 
derstanding of these first clouds is required before def- 
inite conclusions can be made about the formation of 
globular cluster systems. 

In clouds that do survive SN explosions, star for- 
mation halts when there is not enough gas to form 
the shells that trigger the process. Hence, the SN- 
induced star formation method used here is capable 
of self regulation. If these clouds are allowed to ac- 
crete gas from any other source, another sequence of 
star formation is possible. This is consistent with the 
observations that suggest that the Sgr dwarf galaxy 
has had an episodic star formation history (Mighell 
et al. 1999). 

3.2. Stellar Abundances and Metallicities 

The self-enriching fragments of gas, in which the 
first stars are formed, will eventually be destroyed ei- 
ther by internal stellar feedback or by external effects 
such as tidal forces. Star formation is then halted. 
The extremely metal-poor stars left behind contribute 
to the halo field population, and contain informa- 
tion about the nucleosynthesis of elements in zero- 
metallicity stars, their IMF, and the environment in 
which the first clouds formed. Thus, we next com- 
pare the results of our calculation to the observed 
abundance trends in the Galaxy. 

3.3. Metallicities and Ages of Stars 

The metallicities of stars are calculated in this 
model from the metallicities of the shells formed by 
SN events. (In this study, mctallicity always refers 
to the value of [Fe/H].) The first star that explodes 
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is a metal-free Pop-Ill star. Its ejecta depends on its 
initial mass. The ejecta will mix with a given amount 
of primordial gas to determine the metal composi- 
tion of second-generation stars. Figure 2 shows the 
mctallicities produced in shells of different progeni- 
tors as a function of time. They are shown for four 
different progenitor masses: 13, 20, 30 and 40 M©. 
The metallicities produced at early times range from 

[Fe/H] 2.5 to values well below [Fe/H] 4 for the 

most massive progenitors of ~ 40 M©. This implies 
that it might be possible to form Pop-II stars with 
[Fe/H]< —4. This is in agreement with the recent 
discovery of a low-mass star with an iron abundance 
of [Fe/H] = -5.3 (Christlieb et al. 2002). 

Figure 2 also shows that SNell produce shells that 
do not exceed an initial mctallicity of [Fc/H] ^ — 2.5. 
Note, that this uper limit of [Fe/H] < — 2.5 is based 
upon a complete mixing of the supernova ejecta into 
the shell. If this mixing were less efficient, for example 
by the loss of ejecta from the cloud in a jet, this limit- 
ing metallicity would be appropriately reduced. The 
fact that this limit is the value of metallicity at which 
McWilliam et al. (1995) observed a shift in the slope 
for relative abundances vs. mctallicity suggests that 
mixing of ejecta with the shell must be rather efficient. 
This value of mctallicity separates two different stages 
of chemical evolution, the era of inhomogeneous com- 
position dominated by SNell, followed by a stage of 
rapid iron enrichment by type-la SN. As long as the 
mixing of the ejecta is efficient, this also supports the 
idea that only a few SN events are required to raise 
the heavy-element abundances from zero to the values 
observed in the most metal-poor stars in the halo. 

In Figure 3 we show the age-metallicity relation for 
stars in a cloud of initially primordial gas, and assum- 
ing a Salpeter mass function. This cloud was allowed 
to evolve up to 10 Gyr. The ejecta of type-la SN con- 
tributed after a delay time of 1 Gyr. Each point on 
the figure represents a shell that formed stars. The 
mean value of [Fe/H] increases with time, as expected. 
Stars of the same age show a dispersion in metallicities 
that increases for older stars. This is not predicted 
by a simple one-zone model that assumes complete 
mixing, and is a result of the star formation mecha- 
nism applied in this model, i.e. that stars form out of 
the material ejected in Type-II SN events. This will 
allow for a variety of chemical composition depending 
on the progenitor masses and metallicities. 

This large scatter, decreasing with time, can be 
explained by the fact that stars are formed only dur- 



ing SN events. Therefore, stars are the result of the 
mixing of SN ejecta with a limited amount of gas. 
This produces stars of different elemental ratios even 
if they form at the same time, since the SN ejecta are 
dependent on both the initial mass and metallicity of 
the progenitor. 

The most-massive SNe would contribute first and 
reflect their yields in the first Pop-II stars to appear, 
then at later times the products of the less-massive 
SN are expected to appear with different elemental 
ratios (WW95). This does not necessarily mean that 
the first and most massive SNe will produce the lowest 
metallicity Pop-II stars. It is the lowest mass SNell 
that produce shells as low in metallicity as [Fc/H] 
<~ —4, while 20 to 30 M Q progenitors produce higher 
metallicity shells. 

When enough time has passed for the ISM to mix 
on a large scale, the elemental ratios are expected to 
reflect an average of the ejecta of SN events of differ- 
ing masses, including those of type-la. The steep rise 
at early times reflects the dominance of SN individ- 
ual events with different effective yields, and a high 
efficiency for the enrichment of the metal-poor ISM. 
The gradual increase at later times is due to the slower 
enrichment process involving a well-mixed ISM. It is 
clear from the scatter in the figure that a simple age- 
metallicity relation probably does not exist at early 
times. 

3.4. The Metallicity Distribution Function 
(MDF) 

Carney et al. (1996) compared the kinematics and 
chemical properties of metal-poor stars as a function 
of distance from the Galactic plane. They found that 
stars closer to the plane exhibit a slightly more metal- 
rich mean abundance, ([Fe/H]) w —1.7, while those 
farther from the plane have ([Fe/H]} sa —2.0. This 
was postulated to result from the existence of two 
discrete components in the halo. The flattened, inner 
component was probably formed in a manner resem- 
bling the ELS model, while the more spherical, outer 
halo had a large contribution from nearby smaller 
systems similar to the Sagittarius dwarf galaxy or 
other proto-dwarf galaxies. This is in agreement with 
the hierarchical CDM model (Kaufman, White, & 
Guiderdoni 1993). The halo MDF is therefore gen- 
erally believed to be the result of combining the dis- 
tribution functions of individual proto-dwarf galaxies 
accreted by the Galaxy (Cote et al. 2000; Scannapieco 
& Broadhurst 2001, and references therein). 
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Although many recent studies of the formation of 
the halo are supportive of such a scenario (Cote et 
al. 2002; Dincscu 2002; Brook et al. 2003), when 
comparing the observed features of present-day dwarf 
galaxies to those of the (presumably) accreted stars 
and clusters, one finds inconsistencies. The so- called 
satellite catastrophe (Mateo 1998; Moore et al. 1999), 
the angular momentum catastrophe (Navarro & Stein- 
metz 2000), and comparisons of abundances and ages 
in large and dwarf galaxies, lead to the conclusion 
that present day dwarfs are probably not the ma- 
jor building blocks of large galaxies like the Milky 
Way (Shetrone, Cote & Sargent 2001; Baldacci et al. 
2002; Kauffmann et al. 2003; Tolstoy et al. 2003). 
This is not necessarily in disagreement with hierar- 
chical galaxy formation theories. One can still argue 
that the satellite systems that built large galaxies may 
have been different from present-day dwarfs, and that 
these systems are not visible today (Steinmetz 2003). 
For a more detailed discussion, see Tosi (2003). 

The observed MDF of metal-poor stars in the 
Galaxy shows a relatively broad peak with a maxi- 
mum at [Fe/H] ~ —1.6, and a smooth tail extending 
to values of [Fc/H]< — 3. In contrast, the metallicity 
distribution of disk stars shows a more localized peak 
with a sharp cutoff (Ryan & Norris 1991a; 1991b). 
This suggests that the environment in which the disk 
formed was probably more uniform than that which 
formed the halo stars. 

A number of authors (e.g., Norris 1999; Chiba 
& Beers 2000) have suggested that there might be 
a "bump" in the halo MDF at lower metallicities. 
A feature not expected from the simple model, this 
bump might be suggestive of non-uniform enrichment 
at early times. In this picture, pre-Galactic clouds are 
assumed to be the building blocks of the halo. There- 
fore, their metallicity distributions will contribute to 
the total MDF of the halo. The MDF is expected 
to vary from cloud to cloud depending on the initial 
conditions and the duration of chemical enrichment. 

As a crude approximation to the effect of an ex- 
pected range of combined MDF's we made a sim- 
ple study of the MDF's of four different clouds, each 
evolved with a different IMF. We did this under the 
assumption that all these clouds eventually disperse 
and contribute to the field of the halo. The time evo- 
lution of the MDF was followed up to the time of each 
cloud's destruction. In the case of the Salpeter mass 
function, for example, the evolution was followed up 
to 5 Gyr and was found to saturate at 1 Gyr. 



Figure 4 shows the results for the four clouds, 
where the type of mass function chosen for each cloud 
is shown on the figure (all MDFs were normalized to 
unity) . The panels shown in Figure 4 indicate that the 
chemical enrichment in the first cloud, with a simple 
power-law mass function, is much slower than in the 
other three clouds with mass functions biased toward 
massive stars. Therefore the timescales chosen for this 
figure arc different, since the enrichment is slower and 
with the short timescales used in the other three pan- 
els, one would not be able to see the full evolution 
picture in this cloud. 

At the beginning all of the clouds show a disper- 
sion, with the majority of stars at metallicities below 
[Fe/H] = — 5. This describes the first stars that form 
from the shells of the most massive progenitors with 
m > 35 M Q . This could also be concluded from Fig- 
ure 2, which shows that only these massive stars form 
such low- metallicity shells. 

As time goes on, less-massive stars will contribute 
shells in the metallicity range —3 <[Fc/H]< —2.5. 
These produce a rapid increase in the numbers of 
stars in this metallicity range. Moreover, these stars 
should greatly outnumber the first stars formed from 
higher-mass explosions. As soon as stars with m <20 
M start to contribute, shells with metallicities down 
to [Fe/H] ~ —4.5 form and a dispersion is seen once 
again. This dispersion eventually fades out as the 
inhomogeneities start to disappear and the stars be- 
come more metal rich. The final stage shows a peak 
at a different value of metallicity for each cloud. This 
illustrates the different enrichment produced by the 
different mass functions. 

The first cloud, with a simple mass function, achieves 
a peak at [Fe/H] = —2.7, a value consistent with 
the calculations of Ryan et al. (1996) and Nakasato 
& Shigeyama (2000) for the expected metallicity of 
Pop-II stars formed in the shell of a typical SN ex- 
pelled into primordial gas. The other clouds produce 
higher metallicity peaks, corresponding to the larger 
number of massive stars that they produce. These 
peaks have the values [Fe/H]= —2.6, —2.5, and —2.2, 
respectively. The latter corresponds to the fourth 
mass function shown in Figure 1, which produces the 
largest number of massive stars. 

In Figure 5 we illustrate the sum of the final 
stages of all four clouds. This produces a broader 
distribution extending from about [Fe/H]= —3.2 to 
about [Fe/H] = —2.2 and peaked at a value of [Fe/H] 
= —2.6. The MDF observed in the halo shows a yet 
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broader distribution, than that indicated on Figure 
5. This is probably due to the large number of dif- 
ferent systems that contributed to it, and a higher 
value for the maximum mass indicating the contri- 
butions from systems that may have produced stars 
from pre-enriched gas. 

4. Relative Elemental Abundances 

4.1. The Observed Trends in Halo Stars 

Very metal-poor halo stars show a great diversity 
in their absolute elemental abundances and can show 
considerable scatter (up to a factor of 100) in the 
observed abundances of heavy elements relative to 
Fe. (Ryan & Norris 1991; Gratton & Sneden 1994; 
McWilliam et al. 1995; Ryan et al. 1996; McWilliam 
1998; Sneden et al. 1998; Norris, Ryan, & Beers 2001, 
Carretta et al. 2002). At higher metallicities the scat- 
ter gradually decreases until this ratio terminates at 
a value that corresponds to an average over the IMF 
of the element to iron ratios of the stellar yields. 

McWilliam et al. (1995) detected a shift in the 
slope of the abundances of the elements Al, Mn, Co, 
Cr, Sr, and Ba relative to Fe below [Fe/H]~ —2.5. 
They suggested that this unusual chemical composi- 
tion must have originated from the fact that super- 
nova yields must have changed with time, or in other 
words, are metallicity dependent. These results were 
confirmed later by Ryan et al. (1996) who also showed 
that scatter in the abundance ratios increases with 
decreasing [Fe/H]. 

Under the assumption of efficient mixing of the 
supernova ejecta with the clouds, these observations 
support the idea that the most metal-poor stars ex- 
hibit the ejecta of very small numbers of super- 
novae. The metallicity-dependent yields used in our 
present model, together with the unique method of 
SN-induced star formation, are capable of explain- 
ing these trends. The observed relative abundances 
in the halo are fully reproduced with this model for 
several elements. These results are summarized and 
discussed in the following sections. 

4.2. Overview of Model Predictions 

To obtain a picture of the early elemental abun- 
dance evolution, simulations were run for 5 Gyr us- 
ing a massive cloud (10 7 M ), and a Salpeter initial 
mass function. After 1 Gyr, SN of type-la were al- 
lowed to contribute their ejecta. The elemental abun- 



dances were calculated and recorded for every shell 
that forms during the simulation. The model pro- 
duces a small number of stars at very low metallic- 
ities, [Fc/H]< —3, showing a considerable spread in 
[X/Fe] ratios, ranging from less than 0.5 dex in the 
case of Ca and Ti to as much as 1 dex for Mg (Co- 
hen et al. 2004) [see however Arnone et al. (2004) 
who obtain almost no dispersion for Mg]. The ob- 
served scatter of the abundances for model stars is 
given by the spread in metallicities of the SN models. 
Therefore, the large scatter in the abundance ratios 
observed in low-metallicity stars is inherently repro- 
duced by the stochastic nature of this model. 

Local inhomogeneities start to disappear at —3.0 < 
[Fe/H] < —2.5, corresponding to ~ 0.1 Gyr from the 
onset of star formation. At this stage most of the 
massive SNe have contributed, and their ejecta have 
already mixed with the gas in the ISM. This metallic- 
ity range marks the end of the early phase and the be- 
ginning of the transition to the well-mixed phase. At 
this stage, the gas becomes more metal-rich. Hence, 
newly formed stars will no longer exhibit abundance 
patterns of a single SN, but rather an average over 
the IMF of the supernovae that contributed to the 
enrichment of the local ISM. These values resemble 
the predictions of a simple one-zone model in which 
stars are assumed to form out of the well mixed gas 
in the ISM. The spread in the relative abundances de- 
creases gradually, reflecting the ongoing mixing pro- 
cess as more SNe pollute the ISM, and the values of 
[X/Fe] approach the solar value. 

4.3. [a/Fe] 

Figure 6 shows [a/Fe] vs. [Fe/H] for several al- 
pha elements calculated by a straightforward appli- 
cation of the model (small dots). These are com- 
pared with two sets of observational data, Ryan et 
al. (1996) (triangles) and Norris et al. (2001) (circles). 
The observed values show an over-abundance relative 
to solar, and a scatter that is significant at all values 
of [Fc/H], larger at the lowest abundances. There is 
also a peculiar trend in the data to show some stars 
which appear to have increasing [a/Fe] with decreas- 
ing [Fe/H] which shows up as an extension in the data. 
This behaviors is apparent in the abundances of Mg, 
Ca, and Si. Ti, on the other hand, shows a scatter 
that is almost constant over the metallicity range. 

The over-abundance in alpha elements relative to 
iron compared to the Sun can be understood as the 
result of the dominance of SNell at the time, since 
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they are believed to be the major production sites for 
alpha elements. This also indicates that the corre- 
sponding stars formed within the first Gyr of the on- 
set of significant star formation. After 1 Gyr, type-la 
SN events are expected to shift the trends with large 
amounts of Fe. Therefore, the ratios decline gradu- 
ally towards solar metallicities. The so-called plateau 
for [a/Fe] ratios at low metallicities has been shown 
(Arnonc ct al. 2003; Cayrel et al. 2004; Cohen ct 
al. 2004) to have less dispersion (at least for some 
elements) than previously thought (Ryan ct al. 1996; 
McWilliam 1997; Norris et al. 2001). This dispersion 
can set strong limits on the nature of the first SN, and 
the nature of the enrichment process when compared 
to the predictions of chemical evolution models. 

More recent studies of the abundances of metal- 
poor stars show similar trends in general (Carretta et 
al. 2002; Fulbright 2002; Ivans et al. 2002;). In ad- 
dition, all three studies report the detection of stars 
with low alpha-element abundances relative to iron, 
in contrast to the majority of metal-poor stars. To 
explain these observations, they suggest that these 
stars probably witnessed the enrichment of type la SN 
events (Ivans et al. 2002), owe their origins to lower 
mass stellar systems with slower enrichment histories 
(Carretta et al. 2002), or were formed in the Galaxy 
but at times of nonuniform mixing where a mass func- 
tion that produces less high mass stars (which are re- 
sponsible for producing Mg as less massive SNII pro- 
duce the heavier alpha elements) would be respon- 
sible for the low-alpha stars. This is an interesting 
phenomena that needs to be included in a complete 
picture for the chemical evolution of the Galaxy, once 
enough data is available, since it does have the poten- 
tial to set limits on the type and number of SN events 
that took place in the environment from which these 
stars formed. 

Carretta et al. (2002), Arnone et al. (2004), Cayrel 
et al. 2004. and Cohen et al. (2004) have analyzed a 
sample of extremely metal-poor stars, down to low 
metallicities, [Fe/H]> —3.0. They find that the scat- 
ter in the abundance ratios [X/Fe], is surprisingly 
small, even at the lowest metallicities, compared to 
what is expected from a stochastic early evolution. 
They suggest several ideas to reduce the scatter pre- 
dicted by chemical evolution models, including effec- 
tive mixing to maintain the chemical homogeneity in 
the clouds, limiting the mass range of the IMF, intro- 
ducing a first generation of very massive stars (> 100 
M ), and relaxing the closed-box approximation. 



As seen in Figure 6, our most naive model repro- 
duces the basic observational trends for alpha ele- 
ments relative to iron. However, the model tends to 
overestimate the scatter for alpha elements, and ex- 
hibits over all lower values than expected for [X/Fe] 
ratios. This suggests an over production of Fe. This is 
most evident in Mg and Ti, which show mean values 
of [X/Fe] that are 0.5 dex lower than the data all the 
way down to the lowest metallicities. Ca and Si show 
better agreement below [Fe/H] < —3, but also decline 
and reach subsolar values as [Fe/H] increases. The 
suggestion by Timmes et al. (1995) of reducing the 
amount of Fe produced by massive stars seems to be 
required, especially in the cases of elements Mg and 
Ti. This indicates the uncertainty in the Fe yields, 
which is representative of Fe-group elements in gen- 
eral, and will be discussed in the next section. The 
overproduction of Fe can also be attributed to the 
contributions of SNIa. A delay time longer than 1 
Gyr may be required. Although this is a possibility, 
it still seems more likely that the Fe yields from SNell 
are responsible, since the apparent underproduction 
of alpha elements seen in our calculation is not seen 
for Fe group elements as well. 

Figure 7 shows the results of our model when the 
Fe yield is reduced by a factor of 2. The agreement be- 
tween the data and the calculation in this case is much 
improved. We even produce the bifurcation of trends, 
e.g. in [Mg/Fe] and [Si/Fe] there are two trends in the 
data. One is a flattening at [X/Fe]~ 0.5, and another 
a linear increase in [X/Fe] with decreasing metallic- 
ity. In our model, these two trends are an artifact of 
which type of star produced the first supernova. The 
metallicity-dependent yields of Ca reproduce the val- 
ues of [Ca/Fe] and the scatter at low metallicities very 
well, and show a flattening with decreasing metallic- 
ity. This was also observed in the samples analyzed 
by Carretta et al. (2002), Cayrel et al. (2004) and Co- 
hen et al. (2004). Carretta et al. (2002) also suggest 
that the scatter in [Si/Fe] ratios at low metallicity is 
due to uncertainties in the observational techniques 
used for this element. Our results for the alpha ele- 
ment Ti shows lower values than the data and rather 
little dispersion. This is consistent with the most re- 
cent observations (e.g. Cayrel et al. 2004; Cohen et 
al. 2004). However, as for any element, the less abun- 
dant it is, the more uncertain is its production rate, 
as calculated from supernova models. This is surely 
the case tith Ti. 

Previous studies (e.g., Argast et al. 2002) have 
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found that such stochastic models tend to predict too 
much scatter (e.g. for O and Mg) and fail to repro- 
duce the trends of other elements. Our model is also 
based upon a stochastic picture and also can overes- 
timate the scatter in Mg compared to the most re- 
cent observations unless a modification in the yields 
is adopted as in Argast et al. (2002). This will be 
discussed in more detail in §4.5. 

4.4. [Fe-peak / Fe] 

For Fe-peak elements, the element abundance ra- 
tios relative to iron are close to solar down to a 
metallicity [Fe/H]> - 2.5. That is, except for Mn, 
which is underabundant in all halo stars. Below this 
([Fe/H]< - 2.5), the behavior changes (McWilliam et 
al. 1995; Ryan et al. 1996; McWilliam 1997; Norris 
et al. 2001). Both Cr/Fe and Mn/Fe decrease to sub- 
solar values, Co/Fe increases to super-solar values, 
and Ni maintains an average value slightly above so- 
lar. The scatter increases towards lower metallicities. 
The correlation seen between Cr/Fe and Mn/Fe ra- 
tios (both increasing with increasing metallicity), was 
speculated by Carretta et al. (2002) to be an indica- 
tion of the absence of very massive stars (m > 100 
M Q ), since they argue that these stars would under- 
produce nuclei with odd nuclear charge. Therefore, 
[Mn/Fe] would decrease with metallicity while [Cr/Fe] 
would remain approximately constant. They sug- 
gested an explanation for the trends of these ratios 
decreasing with decreasing [Fe/H] as due to variation 
in mass cuts in SNell events as a function of progen- 
itor mass. The observed trends can be reproduced 
if the mass cut is smaller for larger progenitor mass, 
which presumably were more common at lower metal- 
licities. 

Figure 8 shows the calculated values of relative 
abundances for Cr, Co, Mn and Ni. The observed 
ratios of Mn and of Ni are reproduced well by the 
metallicity- dependent yields of WW95 in this model. 
Cr, on the other hand, does not seem to produce 
satisfactory results. The stars produced at very low 
metallicity show a small over-abundance relative to 
solar, contrary to the observed under-abundance. The 
excess in Co/Fe was produced fairly well with the 
metallicity-dependent yields below [Fe/H] = —2.5 if 
the contributions from SN more massive than 40 M Q 
are excluded (see section 4.5). This behavior of Co 
at the lowest metallicities has not been produced by 
any previous models, and no known stellar yields were 
able to predict it even by modifications in explosion 



parameters of SN-II models (e.g. Nakamura et al. 
1999). 

Our calculations raise the questions as to how well 
the nucleosynthesis of massive stars is understood, es- 
pecially for progenitors of low metallicity. The depen- 
dence of the yields of these elements on the position 
of the mass cut and degree of mixing for material 
ejected from SNe is one of the primary sources of large 
errors in calculated abundance yields. Mn and Cr 
are produced mainly during explosive Si burning, and 
therefore have a complicated dependence on the pro- 
genitor mass. Efforts to understand the abundance 
trends of iron-peak elements, which are formed close 
to the mass cut, have focused on the possible alpha- 
rich freeze-out (McWilliam et al. 1995), the location 
of the mass cut, and the dependence of yields on the 
star mass, metallicity, and neutron excess (Nakamura 
et al. 1999). Our model predicts a scatter in the ra- 
tios for Fe-group elements that is larger at very low 
metallicities than the observed values. Mostly this 
large scatter is contributed from the products of SNe 
more massive than 35 M . This may suggest that 
these stars do not contribute and end their lives as a 
black hole (see §4.5). This could also be attributed to 
the several approximations made in the model, espe- 
cially the exclusion of any overlap or mixing among 
SN shells. 

4.5. Contributions from Different Supernovae 

The ejecta of the most massive stars are expected 
to dominate the earliest phases of the chemical evo- 
lution of each cloud. These stars produce values of 
[Fc/H]< -2.5 in their shells. According to WW95, 
at zero metallicity, the most massive stars (35 M Q < 
M < 40M Q ) produce metallicities lower than that 
observed while the lowest observed metallicities of 
[Fe/H] ~ —4 are produced by progenitors of masses 
less than 20 M©. These various behaviors of SN 
with different progenitor masses is responsible for the 
dispersion in metal abundances at low metallicities. 
Studying the dispersion in the relative abundances of 
elements in the halo will allow for constraints to be 
set on chemical evolution models. It can help us un- 
derstand the nature of the first SNe, and the inhomo- 
geneous enrichment processes that took place. It also 
constrains the nature of the IMF and star formation 
rates at the corresponding times. 

Figures 9 to 12 show the contributions to the el- 
emental abundances by different progenitor masses 
down to metallicities of [Fe/H]= —5 and up to —2.5. 
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This range represents the early stages of chemical evo- 
lution in the Galaxy for which our model is applicable. 
Once [Fe/H] > —2.5, the contributions from SNIa pol- 
lute the ISM and the contributions from SN of Type 
II become less dominant. For these figures, we have 
run the simulation allowing only a range of progenitor 
masses to contribute its ejecta. The mass ranges are 
shown on the figures. The same is shown for Fe-group 
elements in Figures 13 to 16. 

The main feature that can be seen in all eight fig- 
ures is that the progenitors responsible for producing 
extremely low metallicities in the Galaxy are either 
more massive than 35 M Q or less massive than 20 Mq . 
Figure 2 shows that these extremely low-metal stars 
are only produced at the earliest times, and therefore, 
are probably the products of Pop-Ill stars. Also, it 
is the the high-mass progenitors m>40 M Q that pro- 
duce shells with metallicities way below [Fe/H] = — 5, 
while the progenitors with m<20 Mq can only pro- 
duce shells with metallicities down to [Fe/H]> — 4. 
Therefore, the fact that the ejecta of the high mass 
progenitors is missing at low metallicities, while the 
current belief is that they have been available in the 
halo from the beginning, can definitely be attributed 
to a lack of data. And while the search for lower 
metallicity stars still continues, this will remain an 
open question. 

If the extremely low-metal stars are never found 
in our Galaxy, one can suggest that the first stars 
that formed in pre-Galactic clouds formed out of pre- 
enriched gas, and that there are no Pop-Ill stars in 
our Galaxy. Or that stars more massive than 35 M Q 
do not contribute their ejecta, and instead, end their 
lives as black holes. 

The missing ejecta of SN events with progenitors 
less massive than 20 M Q in the data, at low metal- 
licities, may be suggestive of the fact that such stars 
are not produced at the earliest times, in support of 
the theory that a top-heavy IMF dominated star for- 
mation in pre-Galactic clouds. Another feature in 
the figures supporting this idea is that for several el- 
ements (Mg, Cr, Co, Ni and Mn), the calculated rel- 
ative abundances that lie outside the observed ranges 
at metallicities [Fe/H] > —4 are either contributed 
partially (as in the cases of Ni, Cr, and Co), or totally 
(as in the cases of Mg and Mn) by stars less massive 
than 25 M . This is suggestive of the fact that SNe 
that contributed before [Fe/H] ~ — 2.5 were all mas- 
sive. In the cases of the elements Ni, Co and Cr, we 
can also see contributions outside the observed ranges 



of relative abundances coming from stars more mas- 
sive than 35 M Q . This also may imply that these stars 
do not contribute as SNell. When these SNe at the 
low-mass and high-mass end are not allowed to con- 
tribute at the earliest stages, the model fits very well 
the recent observational data (Caretta et al. 2002; 
Cayrel et al. 2004; Arnonne et al. 2004) which report 
less scatter than previous studies. 

Again we need to keep in mind here that the nu- 
cleosynthesis calculations of these elements are still 
not certain enough to allow us to build accurate con- 
clusions regarding their production. Therefore, we do 
see here that there is a great need to study these early 
stages of chemical evolution in more detail in order to 
be able to set limits on the details of the production 
of these elements. 

The calculated values from the high-mass progen- 
itors do not fit the observed values for Co even at 
higher metallicities, while the rest of the mass range 
down to 10 M Q reproduces the observations very well. 
This result for the low-mass end is in contrast to what 
is found for the other elements Mg, Ti, Cr, and Mn. 
In these cases, the low-mass progenitors produce val- 
ues of the relative abundances out of the range of the 
data. This may be suggestive of the fact that the 
first stars were massive and that the SN that con- 
tributed before [Fe/H]~ —2.5 were all massive SN. 
The alpha-element enhancement observed for stars 
with [Fe/H]< —2.5 is indicative of a high-mass IMF 
in the first Gyr of the Galaxy's formation. After 
this time, as type-la SN contribute their Fe ejecta, 
at [Fe/H]~ —2.5, the alpha-element ratios go down 
to solar values. 

5. Summary and Conclusion 

We developed a CDM-motivated model to simu- 
late the stochastic early chemical evolution of proto- 
Galactic clouds in the halo. This model is based upon 
a SN-induced star-formation mechanism which keeps 
track of the chemical enrichment and energy input to 
the clouds by Type II and Type la supernovae. An 
important feature of this model is the implementa- 
tion of metallicity-dependent yields for all elements 
at all times, and the inclusion of finite stellar life- 
times. These were found to have important effects 
on the metallicity distribution functions for stars in 
the clouds, their age-metallicity relation, and relative 
elemental abundances for of alpha- and Fe-group ele- 
ments. 
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The stability of these clouds against destruction 
was considered, and results were found to depend 
upon the choice of initial mass function for low- 
mctallicity stars. Most clouds with initial masses sim- 
ilar to presently observed globular clusters were found 
to survive disruption from the onset of star forma- 
tion, suggesting that these systems could have been 
responsible for the formation of the first stars, and 
could provide an environment for the self-enrichment 
of globular clusters. However, such clouds are only 
stable when one assumes an initial mass function that 
is not too biased towards massive stars, indicating 
that even if the first stars were formed according to 
a top-heavy mass function, subsequent star formation 
may need to have proceeded with a present- day mass 
function, or happened in an episodic manner. 

An important conclusion of this study is that the 
dispersion of the metallicity distribution function ob- 
served in the outer halo is naturally reproduced by 
contributions from many clouds with different initial 
conditions. The scatter in metallicity as a function of 
age for stars in our simulations is very large at early 
times, implying that no age-metallicity relation exists 
in the very earliest stages of galaxy formation. 

Regarding observed elemental abundances, the ob- 
served scatter is reproduced fairly well for most, but 
not all, elements. In particular, the predicted relative 
abundances of alpha elements compared to iron can 
be made to agree with the observed values down to a 
metallicity below [Fe/H] ~ —4, but only if the total 
iron produced is reduced by a factor of two from our 
adopted WW95 yields. Similarly, iron group elements 
Cr, Co, Mn, and Ni all exhibit abundance trends vs. 
[Fe/H] that are somewhat different than those ob- 
served. The simplest explanation for both of these 
phenomena within our model is to move the mass-cut 
for iron and iron-group elements outward in the stel- 
lar models for all stars except those with masses in 
the range 25 < m < 35 M©. If yields are restricted in 
this way a good fit to all elements can be obtained. 
We note, however, that though this is a simple expla- 
nation, we can not rule out that these trends could 
result from some of the previously noted shortcomings 
of the model rather than a stellar evolution effect. For 
example, the reduced iron and iron-group yields could 
perhaps equally well have been achieved by relaxing 
the assumption of efficient mixing of the ejecta with 
the cloud. 

A particularly reassuring result of the present study 
with regard to alpha elements is that these models 



naturally predict the observed detailed trends. That 
is, these models reproduce both stars which increase 
in [a/Fe] with decreasing [Fe/H] and stars which ob- 
tain a near constant value of [a/Fe] albeit with a 
significant scatter. In these models, this bifurcation 
of the abundance distributions can be attributed to 
whether the initial supernovae within the clouds were 
massive m ~ 40 M Q progenitors or of lower mass. 

The contributions to the abundances from super- 
novae with different progenitor masses and metallic- 
ity suggest that the low-mass end of Type-II SN was 
probably absent at the very lowest metallicities, and 
that the upper mass limit for the first stars that con- 
tributed to nucleosynthesis may be <40 M©. 
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Fig. 2. — Metallicity enrichment vs. time in the first 
2 Gyr for shells containing progenitor stars of various 
masses. The progenitors shown have masses 13, 20, 
30 and 40 M Q 

Fig. 3. — Age-metallicity relation for stars produced 
by the model. Note the large dispersion at early 
times. 

Fig. 4. — The metallicity distribution functions pro- 
duced in clouds with four different initial mass func- 
tions as labeled. 

Fig. 5. — The sum of the MDFs for the four clouds 
shown in figure 4 

Fig. 6. — The abundances of the alpha elements Mg, 
Ca, Si, and Ti. Data is from Ryan et al. (1996) (tri- 
angles) and Norris et al. (2001) (circles). The small 
(dots) represent the model stars 

Fig. 7. — The abundances of the alpha elements when 
Fe yields are reduced by a factor of 2. Data is from 
Ryan et al. (1996) (triangles) and Norris et al. (2001) 
(circles). The small (dots) represent the model stars 

Fig. 8. — The abundances of the Fe-group elements, 
Mn, Ni, Cr, and Co . Data is from Ryan et al. (1996) 
(triangles) and Norris et al. (2001) (circles). The 
small (dots) represent the model stars 



Fig. 9. — The contribution from different progenitor 
masses to the value of [Ca/Fe]. The mass ranges are 
given in solar masses. The top left figure shows the 
whole range. Data is from Ryan et al. (1996) (tri- 
angles) and Norris et al. (2001) (circles). The small 
(dots) represent the model stars. 
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Fig. 10. — The contribution from different progenitor 
masses to the value of [Mg/Fe]. The mass ranges are 
given in solar masses. The top left figure shows the 
whole range. Data is from Ryan et al. (1996) (tri- 
angles) and Norris et al. (2001) (circles). The small 
(dots) represent the model stars 
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Fig. 11. — The contribution from different progenitor 
masses to the value of [Si/Fe]. The mass ranges are 
given in solar masses. The top left figure shows the 
whole range. Data is from Ryan et al. (1996) (tri- 
angles) and Norris et al. (2001) (circles). The small 
(dots) represent the model stars 



Fig. 12. — The contribution from different progenitor 
masses to the value of [Ti/Fe]. The mass ranges are 
given in solar masses. The top left figure shows the 
whole range. Data is from Ryan et al. (1996) (tri- 
angles) and Norris et al. (2001) (circles). The small 
(dots) represent the model stars 



Fig. 13. — The contribution from different progenitor 
masses to the value of [Cr/Fe]. The mass ranges are 
given in solar masses. The top left figure shows the 
whole range. Data is from Ryan et al. (1996) (tri- 
angles) and Norris et al. (2001) (circles). The small 
(dots) represent the model stars. 



Fig. 14. — The contribution from different progenitor 
masses to the value of [Co/Fe]. The mass ranges are 
given in solar masses. The top left figure shows the 
whole range. Data is from Ryan et al. (1996) (tri- 
angles) and Norris et al. (2001) (circles). The small 
(dots) represent the model stars 



Fig. 15. — The contribution from different progenitor 
masses to the value of [Mn/Fe]. The mass ranges are 
given in solar masses. The top left figure shows the 
whole range. Data is from Ryan et al. (1996) (tri- 
angles) and Norris et al. (2001) (circles). The small 
(dots) represent the model stars. 



Fig. 16. — The contribution from different progenitor 
masses to the value of [Ni/Fe]. The mass ranges are 
given in solar masses. The top left figure shows the 
whole range. Data is from Ryan et al. (1996) (tri- 
angles) and Norris et al. (2001) (circles). The small 
(dots) represent the model stars 
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